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ABSTRACT 

Magnetic helicity has risen to be a major player in dynamo theory, with the helicity of the small-scale 
field being linked to the dynamo saturation process for the large-scale field. It is a nearly conserved quantity, 
which allows its evolution equation to be written in terms of production and flux terms. The flux term can be 
decomposed in a variety of fashions. One particular contribution that has been expected to play a significant role 
in dynamos in the presence of mean shear was isolated by Vishniac & Cho (2001, ApJ 550, 752). Magnetic 
helicity fluxes are explicitly gauge dependent however, and the correlations that have come to be called the 
Vishniac-Cho flux were determined in the Coulomb gauge, which turns out to be fraught with complications 
in shearing systems. While the fluxes of small-scale helicity are explicitly gauge dependent, their divergences 
can be gauge independent. We use this property to investigate magnetic helicity fluxes of small-scale field 
through direct numerical simulations in a shearing-box system and find that in a numerically usable gauge the 
divergence of the small-scale helicity flux vanishes, while the divergence of the Vishniac-Cho flux remains 
finite. We attribute this seeming contradiction to the existence of horizontal fluxes of small-scale magnetic 
helicity with finite divergences. 

Subject headings: MHD — turbulence — Sun: magnetic fields 



1. INTRODUCTION 

The large-scale magnetic field of the Sun and 
other stars is often mode led using mean-field theory 
jRiidiger & Hollerbachll2004h . Important ingredients in this 
theory are the a effect responsible for field amplification and 
an enhanced (turb ulent) magnetic diffusivity MofFatt (1978j); 
iKrause & Radlerl (1198 0). Once the field has reached apprecia- 
ble field strength, these effects become modified through the 
backreaction of the Lorentz force. Often a simple algebraic 
quenching formula is being assumed, but such a simple 
prescription is unable to model corr ectly the quenching under 
more general conditions with shear (Brandenburg et all200Ih 
or boundaries (Brandenburg & Dobler 2001). 

Over the past decade, significant progress has been made in 
modeling the dynamo saturation process in mean-field mod- 
els through the development and use of the dynamical a 
quenching methodology. This methodology (originally due 
to Kleeorin & Ruzmaikin 1982) has explained several puz- 
zling features of MHD dynamos, such as the slow satura- 
tion phase of a homogeneous dynamo (Field & Blackmaiil 
120021: iBlackman & Brandenburg! I2002h . and has reposed 
the crucial question of catastrop hic a quenching (see 
[Brandenburg & SubramanianI 12005 a l for a review). In this 
picture, the saturation of the dynamo is caused by the build- 
up of magnetic helicity, which is nearly conserved in the high 
conductivity limit in the absence of fluxes. This raises the 
possibility of speeding up the saturation process and reaching 
significant saturation field strength through mechanisms that 
export or destroy magnetic helicity. 

Making general use of the dynamical a quenching metho- 
dology in open systems then requires an understanding of 
magnetic helicity fluxes, more significantly an understand- 
ing of the fluxes of magnetic helicity of the small-scale field. 
In the following we often refer to such fluxes as small-scale 
magnetic helicity fluxes, although this is not quite accurate, 

|alex.i.hubbard@gmail.com (Revision: 1.81 J] 



because it is itself a mean quantity and not a fluctuation. Re- 
cent work has shown that in inhomogeneous systems there 
is a turbulent diffusive flux of small-scale helicity, at least 
in the absence of shear, although the diff usion coefficient 
is in some cases smaller than expected ( Mitra et al.l 120 lOt 
iHubbard & Brandenburg! I2OIOI) . Consideration of that flux 
term has allowed mean-field models to capture the saturation 
behavior of some non-triply periodic, non-homogeneous dy- 
namos without shear Unfortunately, the fluxes of small-scale 
magnetic helicity are explicitly gauge dependent, and in the 
presence of turbulence, can be decomposed in diflrerent fash- 
ions. ^ 

In IVishniac & Choi (12001!) . an interesting component of 
the flux of small-scale helicity was isolated. This com- 
ponent has been named the Vishniac-Cho flux, henceforth 
the yC flux. In later w ork ( Subramanian & Brandenburp 
120041: [Brandenburg & Subramanian 200 5b). the form of this 
flux in the presence of uniform shear was calculated, and 
found to be both simple and of significant magnitude. Shear 
drives an D. effect and is an important and nearly omni- 
present player in astrophysical dynamos, so the VC flux 
has seen significant interest, both in mean-field modeling 
dBrandenburg & SubramanimI l2005bb and in the interpreta- 
tion of the diff'erences between similar direct n umerical sim- 
ulations with diff^e ring boundary conditions dBrandenburg! 
12005!: iKkpyla et an i2C)08). It is important to realize however 
that shear poses unique difficulties in the formulation, and 
importantly, interpretation, of magnetic helicity fluxes. It is 
the goal of this paper to explore those difficulties and deter- 
mine the significance of the VC flux. We use the perhaps 
surprising result that while magnetic helicity fluxes are gauge 
depende nt, their divergen ces may not be across broad gauge- 
families (M itra et al.l2 010) to allow us to compare the VC flux 
with the small-scale magnetic helicity flux in a gauge where 
the mean shear is easy to treat. Our investigations will bear 
weight on the interpretation and use of the VC flux, but we 
will not and indeed cannot extract the VC flux from the simu- 
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lations we perform. 

In Section |2] we discuss the various contributions to the 
magnetic heHcity flux and define our mean-field decomposi- 
tion. Further, we explain the broad gauge independence of 
small-scale magnetic helicity flux divergences. In Section [3] 
we sketch the difficulties inherent in uniform shear, define 
the shearing-advective gauge and derive the small-scale mag- 
netic helicity flux in that gauge. In Section |4] we discuss the 
VC flux as calculated injSu bramanian & Brandenburg ( 2004) 
and iBrandenburg & Subramaniani C2005b). In Section \5\ we 
present the results of our direct numerical simulations and 
compare the results with the VC flux. We discuss the signifi- 
cance of our results in Section|6]and conclude in Section|2l 

2. MAGNETIC HELICITY FLUX AND FORMALISM 

We begin by deriving the formula for magnetic helicity 
fluxes in general. The MHD equations for the magnetic field 
are: 



B = Vx A, 
/ = V X B/fio, 
E^-UxB + T]fHiJ, 
dA 

— = - VO, 

dt 



(1) 
(2) 
(3) 

(4) 



where rj is the molecular resistivity, /io is the vacuum per- 
meability, and <E) is the electrostatic or scalar potential that 
determines our gauge. For example, setting O = results 
in the Weyl gauge, of interest numerically because it sim- 
plifies Equation (|4]i, while a solution of Equation (|4]i with 
V^O = -V S will have constant V- A and with an appropriate 
initial condition on A results in the Coulomb gauge. 

The time evolution of the magnetic hehcity density h = A-B 
is then given by 



dh dA 

B 

dt dt 



■ A ■ -— 
dt 



= -2tjhoJ B -V (Ex A + <DB) . 



(5) 



The flux f o of magnetic helicity in a given gauge with a cor- 
responding <i> can be read out of Equation Q: 

f o = -E X A -kDB = - (17 X B) X A + i&B + //yUo/ x A 
^hU + (<:)-U ■ A)B + T]iJoJ xA. (6) 

In Equation ^ we recognize the advective flux f adv = hU,a 
resistive flux f res = Wo J x A, and finally a dynamical flux, 

f dyn = (CD - ■ A)B. (7) 

The formula for f dyn leads us to consider "advective" gauges 
of the form O = t/' ■ A, where U' is one component of the 
velocity field, see Section [STI 

2.1. Mean-field decomposition 

We proceed to a mean-field decomposition of the magnetic 
helicity flux. We denote general averaging schemes by over- 
bars. Fluctuating terms will be denoted by lower cases or 
primes: 



A = A -H a, A B ^A B + (A B)' . (8) 
The mean-field decomposition of Equations (|4| and (|5]l yields 



dA - - - - - 
— ^UxB + G- TjfioJ - VO, 

dt 



where E> = uxb and 
dh 



— = -2riLH)J B-V ExA-V Q>B. (10) 
dt 



The latter can be written as 



dh dhm dhf 
'dt ~ ~dt' 'di' 



(11) 



where h - h^ + hf, with h^ = A ■ B being the helicity in the 

large-scale fields and hf = a ■ b the helicity in the small-scale 
fields. 

The evolution equations for these helicities are 



+ 2G B- 2rim)J -B-V ■[ExA + (DB) , (12) 



dhm 
dt 

^ ^-2GB-2i]HQPb-V ■(el<a+'^], (13) 



dt 

where = (D - <1) is the fluctuating scalar potential and 
£>- aB - j]tixoJ 



(14) 



is the mean turbulent electromotive force, where a is the a 
effect and % is the turbulent magnetic diffusivity. From Equa- 
tions ( fT2l i and ( fTsT l we find the fluxes of the large-scale and 
small-scale fields: 



Fm^Ex A + OB, 
F[-e xa + tpb. 



(15) 
(16) 



(9) 



2.2. The significance of gauges for magnetic helicity fiuxes 

It is clear from the form of f dyn in Equation (|7| that any 
consideration of helicity fluxes must also take into account 
the gauge choice, as at any point where Fdyn 0, the value 
of f dyn can be arbitrarily set by the gauge or, equivalently, the 
condition for Fdy„ being independent of O is that f dyn = 
for all gauges. A gauge choice that generates a desired flux 
along X is always possible provided B ■ x Q. Such a gauge 
takes the form <l)(jc, f) = fix, t)B^. Recall also that boundary 
or symmetry conditions on the physical system do not apply 
to the vector potential or the gauge (although numerical simu- 
lations may require gauge choices where they do). The ability 
to add an arbitrary flux of magnetic helicity to the system via 
the addition of a new gauge makes the isolation of differing 
components of the flux a risky business. 

There are effects that mitigate this gauge dependence how- 
ever The divergence of Ff and the term dhjdt are the only 
gauge-dependent terms in Equation ( fT3] ). If hf is indeed gauge 
independent then the divergence of Ff must be gauge inde- 
pendent as well, even though the ffux itself is explicit gauge 
dependent. Clearly the divergence of Ff is the same for all 
gauges for which djif is the same. As long as our shearing- 
box has a time-constant hf in the saturated regime, we can 
make statements about the divergence of Ff for all gauges 
which would have a time-constant saturated /if. 

The dynamical a quenching methodology, one of the pri- 
mary consumers of magnetic helicity information, assumes 
that the small-scale magnetic helicity /if and the small-scale 
current helicity j ■ b are proportional. This requirement is 
often used as an argument in favor of the Coulomb gauge 
dKIeeorin & Rogachevskiiiri999 ). Even if the saturated cur- 
rent helicity is not time-independent (as might be the case 
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for an oscillating solution), as the current helicity is gauge 
independent, the validity of the dynamic a quenching me- 
thodology assumes therefore that the form er is as well. Re- 
cent studies have supported this hypothesis (iMitra et al.ll2010l : 
iHubbard & Brandenburgll2010 ). at least in the limits of nu- 
merical simulations, which disallow extreme levels of gauge 
pathology by forcing the vector potential to be numerically re- 
solved. (This will be discussed in more detail in a separate pa- 
per where we solve an evolution equation for the gauge trans- 
formation.) Alternatively, one could restrict oneself to fami- 
lies of gauges where the relation holds. We note that, in that 
regard, both the magnetic and current helicities are shearing- 
periodic in our system in the shearing-advective gauge that 
will be described later, while in the Weyl gauge the magnetic 
helicity is not . 

3. SHEAR 

The presence of shear poses further difficulties when con- 
sidering magnetic helicity fluxes. To see this, consider a 
shearing periodic box with an imposed flow. Us - (0,8 x, 0), 
and sides of length L centered on the origin. We will use the 
Weyl gauge (O - 0) and a planar averaging scheme: 

_ pL/2 pL/2 

B(x,y,z,t) = L-^ dxdy'B{x',y',z,t). (17) 

J-L/2 J-t/2 

In what follows, we will assume that the shear flow is the 
only large-scale velocity, and note that it averages to and 
so technically is a fluctuating field. This could in principle 
be avoided under a local planar av erage ov er a square cen- 
tered on {x,y) for example (Brandenburg et al. 2008). How- 
ever, such an average is problematic too, because it does not 
obey one of the Reynolds rules: the average of a product of an 
average and a fluctuation does not vanish. Nevertheless, even 
though uniform shear can complicate averaging schemes, it is 
easier to treat than non-uniform shear with the resulting non- 
uniform Q effect. We therefore proceed with our standard 
(non-sliding) averaging scheme. We will also assume that, 
at the single instant in time that we consider, the vector po- 
tential A is shearing periodic as well. Note that in the Weyl 
gauge, the vector potential will not remain shearing periodic. 
Numerical simulations in the shearing box approximation use 
therefore a different gauge (iBrandenburg et al.il 19951) . as will 
be discussed below. 

The difficulty in treating the helicity flux can be seen from 
Equation (|6]l which, in the Weyl gauge, becomes 

F ^hUs- (Us ■ A)B + hu-(u- A)B + riiiaJ x A. (18) 

If A were shearing periodic, the last four terms of Equa- 
tion ( fTSl l would be likewise shearing periodic and hence 
would not contribute a net divergence to the system. However, 
the first two terms on the right-hand side of Equation (fTsT l vio- 
late shearing-periodicity. In particular, the x component of the 
second term on the right-hand side of Equation (fTsT l reduces 
to 

= -{Us ■ A)B, + ... = -SxAyB, + (19) 

where the dropped terms cannot contribute a net divergence. 
Recall that, at first, AyB^ would here still be shearing peri- 
odic. Systems with shearing-periodic magnetic vector poten- 
tials would then allow for a horizontal flux of magnetic helic- 
ity. 

3.1. Helicity fluxes in advective gauges 



To examine the VC flux numerically, we adopt a homoge- 
neous shearing-periodic setup. (The resulting magnetic field 
will however become inhomogeneous and could produce fi- 
nite magnetic helicity fluxes and flux divergences.) As dis- 
cussed above, to keep the magnetic vector potential itself 
shearing-periodic we must use an appropriate gauge, namely 
= Us ■ A, which we term "shearing-advecti ve". This is 
also the gauge used bv IBrandenburg et alj (Il995h . More gen- 
erally, we can define a family of "advective gauges" with 
- Ua-A, for a component of the velocity Ua, with the cor- 
responding Una - U - Ua- The name "advective" is chosen 
because in this gauge the effect of Ua on the helicity flux is 
advective as can be seen from Equation which becomes: 

Ea^-UnaxB + woJ, (20) 
F^^hUA + (EAxA). (21) 

If Ua is a mean flow (Ua = Ua), then the mean flux of the 
small-scale helicity becomes 

Ff = hfUA + ca xa. (22) 

Alternatively, if Ua is not a mean flow {Ua - 0), then we have 

Ff = h\U'^ + ba xa. (23) 

For our system, with Ua - Us, Sa, a and b are all shearing- 
periodic, and so their mean values, as well as all other mean 
values, are functions of z alone. Correspondingly, only the z 
component of Ff can have a finite divergence, and we have 
eliminated the worry of horizontal magnetic helicity fluxes. 
Further, to the extent that the system is homogeneous, and 
invariant under a 180 degree rotation about the z axis, the hor- 
izontal fluxes vanish entirely except for the advective flux due 
to the shear flow. 

4. THE VISHNIAC-CHO FLUX WITH MEAN SHEAR 

The VC flux, Fvc, has been calculated in several places 
using the first order smoothing approximation and later the 
r approximation. Their applicability to highly turbulent sys- 
tems with large magnetic Reynolds numbers cannot be guar- 
anteed and is subject to verification by numerical simula- 
tions, although it should work in cases of small magnetic 
Reynolds numbers considered in this paper This flux was 
originally calculated in the Coulomb gauge, but it can be 
calculated in a related gauge in which the magnetic helic- 
ity density corresponds to a densit y of magnetic linkages 
dSubramanian & Brandenburg 120061) . It is most interest- 
ing in the case of shear, and we will restrict ourselves to 
the consideration of uniform shear which, as noted above, 
raises concerns about horizontal fluxes with finite diver- 
ge nce. In this s ystem, F yc wa s calculated in the appendix 
of lBrandenburg & SubramanianI (l2005b) to be 

^'vc = Cvc^(fi'-fi'), (24) 

where Cyc is a coefficient expected to be of order unity and h 
is the waven umbe r of the energy-carrying eddies. As eluded 
to in Section 12.21 we assume here that the current helicity is 
proportional to times the magnetic helicity density. There 
are other components of the magnetic helicity flux known, 
for example the term in Equation (fT9] l can be found in the 
F''""' of 'Subra manian & Brandenburd (120061) . and, as noted 
in Section 12.21 the interesting value is the divergence of the 
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total of all fluxes, which might pick up contributions from the 
X component of the flux. 

While the VC flux is of pressing interest to mean-field dy- 
namo theory and it has been invoked in th e interpretation of 
numerical simulations (' Kapvla et al.l I2008h . the work on its 
implications in mean-field theory is not well developed. In 
Section 15.21 we present a brief mean-field analysis of the ef- 
fects of a VC flux in our shearing-sh eet system using (l24ll. a nd 
compare those results to the work of iGuerrero et al.l (120101) in 
spherical shells. 

5. MODEL CALCULATIONS 
5.L Preliminary considerations 

In order to quantify shear-driven magnetic helicity 
fluxes, we consider the shearing box approximation 
(IWisdom & Tremaind 119881) with periodic boundary condi- 
tions in the y and z directions and shearing-periodic boundary 
conditions in the x direction. According to Equation ( l24l l we 
expect a magnetic helicity flux in the z direction. However, 
because our system is periodic in the z direction, there will be 
no net magnetic helicity flux in or out of the domain. Never- 
theless, the local divergence of Fvc should be finite because, 

, —2 —2 

contrary to homogeneous a dynamos without shear, B^ - B,, 
is in general z-dependent for o-Q dynamos. Indeed, the mean 
field that develops is a reasonable approximation to an aO. 
dynamo, where \S \ > \ak-.\. The marginally excited kinematic 
solu tion of a mean-field aO. dynam o is a traveling wave (see, 
e.g.. lBrandenburg & Sokololl2002h . i.e. 

B = Bo (sin 6*, V2 ^ sinC^ + x), o) (25) 



equations 



with 



c - + 



aS 



2k, 



1/2 



' - k~iz - ct), X - "'"4^' 



(26) 



where the sign in front of c is given by the sign of the product 
aS and Bq is an undetermined amplitude factor Note that the 
magnetic helicity of this large-scale field, - \cla\K^B'^, is 
independent of z for the "natural" vector potential 



(-V2 



cos(6' + ;i'), cos 0, 



(27) 



The point of this discussion is to emphasize that even for an 
initially homogeneous system. Equation ^ would predict the 
appearance of a magnetic helicity flux. This flux would lead 
to the annihilation of magnetic helicity fluctuations of oppo- 
site sign — even if such fluctuations were not present initially. 
The effect of such fluxes can be predicted by mean-field mod- 
els w ith catastrophic quenching included (iBrandenburg et all 
l2009h . As we will show in the next section, the VC flux has 
actually an adverse effect on the saturation behavior, and that 
only fluxes with the opposite sign are able to accelerate the 
saturation of the mean field. 

5.2. Mean-field model with difi^usive and VC fluxes 

To demonstrate the difference between mean-field predic- 
tions with and without the presence of a VC flux that is not 
compensated for by other fluxes we use a mean-field dy- 
namical a quenching mo del ( Kleeorin & Ruzmaikin 1982; 
iBlackman & Bra ndenburg! I20b2t iBrand enburg e t al.i i2009l) . 
This methodology combines the dynamical a quenching 



2 

a(z, t) - aK + aM, am - rj^k^ 



daM ^ ,2 



G ■ B um 



Re. 



(28) 



(29) 



with the standard mean-field equation (in the shearing- 
advective gauge with U = 0): 



— - -SAyX + S- 
ot 



WoJ, 



(30) 



where the £ is given by Equation ( fT4b . 
We test diffusive and VC fluxes, setting 

Bf ^ ' 



(31) 



eq 



We measure the strength of the kinetic a effect and the 
shear S with the dynamo numbers 



^ _ oiK ki _ S 

Cff = — — - — , Cs = —J, 
TJiki ki Tjxk-. 



(32) 



where ki - InIL-. is the minimal wavenumber of the domain 
in the z direction. Note that aic is assumed independent of B, 
so the kinetic a effect is therefore also a kinematic one. 

The results are shown in Figure [1] where we plot the satu- 
ration behavior of models with rjjri = 10^, and dynamo num- 
bers Ca = -0.2 and Cs - -20. The left panel covers four 
values of Ka/rji the turbulent diffusion coefficient for hf, with 
the VC flux turned off (Cyc = 0). The right panel displays the 
dynamo behavior for three values of Cyc with Ka/rjt = 0.2. 
Note that in all calculations an early intermediate saturation 
level of ^ 0.15Beq is reached. This is followed by a resis- 
tively slow saturation phase, as was expected from models 
without shear (Brandenburg 2001; Blackman & Brandenburg 
I2OO2I) . In agreement with earlier wo rk, the saturation be havior 
is accelerated by diffusive fluxes (iBrandenburg et af] 12009) . 
In the absence of a diffusive flux, Ka = 0, the field drops sud- 
denly back to lower values and continues to oscillate. These 
oscillations are eliminated by small values of /c„ while not 
significantly affecting the intermediate saturation behavior if 
Ka = 0.2?7t. Note that for positive values of Cvc> the VC flux 
actually has an adverse effect on the saturation behavior and 
only negative values are able to accelerate the saturation. This 
is simi lar to results for mean -field dynamo action in spherical 
shells (IGuerrero et al.ll2010l) . 

5.3. Simulations of shear flow turbulence 

We turn to the computation of magnetic helicity fluxes 
through direct numerical simulations. We solve the stochas- 
tically forced isothermal hydromagnetic equations in a gener- 
ally cubical domain of size {27{f in the presence of a uniform 
shear flow. Us - (Q,Sx, 0), with S - const. 



DA , 
Dt ■ ' 



(33) 



^ = -5 -clVlnp+-JxB + Fvisc + /, (34) 
Ut p 



Dlnp 
Dt 



= -V I/, 



(35) 
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lei machines. Our model set up is identical to that used by 
iKapvla & Brandenburd (l2009h . who studied the frequency of 
dynamo waves in the saturated regime as a function of the 
fractional helicity and thereby the effective dynamo number 
and, more importantly, different magnetic field strengths. 

5.4. The VC flux in simulations 

We focus here on the results of four simulations, with pa- 
rameters given in Table [1] Standard estimates suggest that 
the two dynamo parameters given in Equation ( [32] ) suffice for 
large-scale dynamo action {CaCs > 2), and also their ratio, 
Cs ICa ~ 3Sh kf/ki, is large enough for oscillatory dynamo 
action with significant dynamo wave speeds - even for Run D, 
with a ratio of unity (note the finite wavespeed c - 0.24?7,o^,). 
Recall that the vertical elongation of Run D suppresses non- 
vertical mean field structures by increasing the minimum hor- 
izontal wavenumber to 4, suppressing any potential .jc- varying 
field. 

The magnetic field is normalized to the equipartition value 

fieq = (/^0P0)'^^Mrms, (37) 

while, as suggested by Equation {2M and Section ISTl /if and 
Ff are normalized to 

ho = kf{(Bl}(Bl)f'\ Fo=kfS{W}. (38) 

The brackets represent full volume averaging, here of already 
planar averaged values. 

The simulations developed the expected dynamo wave, so 
we analyze the data in a comoving frame in which the wave 
is standing, allo wing us to average in time. In agreement 
with e arlier work (iMitra et al. 2010; Hub bard & BrandenburgI 

12010), the magnetic helicity of the small-scale magnetic field 
is then statistically steady and therefore the divergence of the 
magnetic helicity flux must be independent of the gauge cho- 
sen. Note that for Run A (Figure|2]i the flux is generally quite 
small (less than lO^^* times the value expected based on equa- 
tion |24] |. except when By = where it shows a small peak. 
The situation is different in Runs B and C (Figures [3] and 

m, where the flux is larger but uncorrelated compared to the 

—2 —2 

spatial dependence expected to be dominated by B,, » B^. 
Finally note that for a mean-field wavenumber k - ki the 
expected VC flux would be of wavenumber 2ki, not much 
smaller than forcing wavenumbers of 3>ki or 5A;i. In Figure|5] 
we therefore present the results of Run D with kf - 2Qki in a 
non-cubic domain with L, - AL. The results are quaUtatively 
similar to Run A which has a similar Re^. 

The significant result to draw from the figures is that the 
small-scale flux is both smaller than the expected Fyc and 
uncorrelated with it. Note that as the y-directed field is much 
greater than the x-directed field, as expected, and that the VC 

TABLE 1 

Summary of the runs discussed in this pa per. The wave speed c is described 
m Equation (26). 

Run Fig. Re^/ Re Sh Ca = k{/ki Cs Resolution c/i^rok. 
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Fig. 1. — Left: Saturation behavior of models with r]t/r] = lO', Cvc = 0, 
Ca = -0.2, and Cs = -20, and for Ka/rji ranging from to 1. Right: Same, 
but for Cvc ranging from -0.2 to +0.2 and Ka/rji = 0.2. 

where D/Dt - d/dt + (U + Us) ■ V is the advective derivative 
with respect to the total flow velocity that also includes the 
shear flow, Cs = const is the isothermal sound speed, f vise - 
p"' V-(2pvS) is the viscous force, = j{Uij+Ujj)-^dijV-U 
is the traceless rate-of-strain tensor, commas denote partial 
differentiation, and / is the forcing term. As in earlier work 
(lBrandenburgl200 ih the forcing function consists of plane po- 
larized waves whose direction and phase change randomly 
from one time step to the next. The modulus of its wavevec- 
tors is taken from a band of wavenumbers around a given av- 
erage wavenumber that is referred to as kf. 

The main control parameters in our simulations are the 
magnetic Reynolds and Prandtl numbers, as well as the shear 
parameter, 

RCM-^, PrM = -, Sh=^. (36) 

We adopt periodic boundary conditions in the y and z direc- 
tions and shearing-periodic boundary conditions in the x di- 
rection. Our initial velocity, in addition to t/s, is t/ = 
and the initial density is p = po = const, while for the mag- 
netic field we take a Beltrami field of negative magnetic he- 
licity and low amplitude (10^^ times the equipartition value). 
The magnetic field grows then exponentially owing to dy- 
namo action and saturates when the field reaches a certain 
multiple of the equipartition value. We solve the governing 
equations using the Pencil Cod^ which is a high-order finite- 
difference code (sixth order in space and third order in time) 
for solving partial differential equations on massively paral- 

' http : //pencil-code . googlecode . com/ 
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Fig. 2. — Panel 1 : visualization of /if as a function of normalized / and z-ct 
for Run A: k[/ki = 5, Re = 0.4, Re^ = 9, and Sli = 0.95. Here c = 0.6f;,o/ti 
is the actual speed of the dynamo wave. Panels 2—4 give the z dependence of 
/if, B, for i = x, y, and F, averaged in the comoving frame over the time of 
the first panel. 

flux is approximately proportional to the square of the By. 
Thus, according to our present results we must conclude that 
there is no support for the validity of Equation (|24] |. 

6. COMPARISON WITH PREVIOUS WORK 

We would like to emphasize that our results are not in 
contradiction with previous calculations: we are not work- 
ing in a gauge where one would expect the VC flux to ex- 
ist. Disentangling the diff ering components (there are four in 
ISubramanian & Brandenb urg 2006, including an unexplored 
triple correlator) is not straightforward. In this work we are 
using the gauge independence of the divergence of the flux of 
small-scale magnetic helicity, as described in Section IZ21 to 
relate the observations in the shearing-advective gauge to the 
expected divergence of f vc ■ 

To date, the only numerical evidence for the VC flux 
comes from interpretations of the differing dynamo behavior 



Fig. 3.— Same as Figure[3] but for Run B: kf/ki = 3, Re = 9, Re^ = 90, 
and Sh = 0.5. In this case, c = 1.047;to^i. The vertical bars in the first panel 
denote the time interval over which the functions in the other 3 panels are 
averaged. Note also that after Hnns^^ff ~ 7800 the dominant dynamo mode 
changes and the field becomes x dependent. 



in shearing systems with vertical field boundary conditions 
(that allow a flux) as compared with those systems with per- 
fect conductor boundary conditions that disallow a flux; see 
Equation Examples o f such indi rect evidence includ e 
the papers by ' Brandenburg (l2005h and'Kap vla et all (120081) . 
An alternative interpretation might simply be that the excita- 
tion condition for the onset of large-scale dynamo action are 
simply delayed sufficiently when changing the boundary con- 
dition from a vertical fie ld to a perfec t cond uctor condition, 
as was discussed also by iKapvla et al.l (120101) . It should also 
be noted that the use of f vc in various dynamical quench- 
ing models has not alleviated catastrophic quenching unless 
Cvc is increased beyond a certain limit where the flux diver- 
gence leads to a magnetic a eflect that is more important than 
the kinematic a effe ct tBrandenburg & Subramanianll2005bl : 
iGuerrero et alJl2OI0l) : see also Appendix]] 
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Fig. 4.— Same as Figure|3] but for Run C: kf jk^ = 3, Re = 19, Re^ = 280, 
and Sh = 0.4. The vertical bars in the first panel denote the time interval over 
which the functions in the other 3 panels are averaged. 

7. CONCLUSIONS 

We conclude that there is at present no evidence for a shear- 
driven vertical flux of small-scale magnetic helicity in a gauge 
where the only significant flux must be vertical. We speculate 
that the finite divergence of the VC flux found earlier in an- 
alytic studies using Coulomb and related gauges might be a 
consequence of the gauge choice, which can generate unex- 
pected horizontal helicity fluxes that are not normally consid- 
ered. When the gauge choice is such that those horizontal 
fluxes are transformed out, there is no remaining vertical flux. 
It appears therefore that the VC flux does either not operate, 
or it at least does not follow the expected functional form. We 
note that diffusive fluxes have been found to exist, so there 
do remain mechanisms that can export small-scale magnetic 
helicity from a dynamo. 

The simultaneous export of large- and small- 
scale helicity at some relative level is inevitable 
felackman & Brandenburg. .2003.) and in fact necessary. 
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Fig. 5.— Same as Figure|3] but for Run D: kfjki = 20, Re = 9, ReM = 7.9, 
and Sh = 0.016. The domain for this run is njl X ;r/2 X 2n, and the hehcity 
in the top panel has been multiplied by 10 to achieve dynamic color range. 
The vertical bars in the first panel denote the time interval over which the 
functions in the other 3 panels are averaged. 

because otherwis e simple considerations ("Brandenbur g et al] 
|200a [Brandenb urg & Subr amanian 2005a) have long sug- 
gested that the magnetic energy would reach unrealistically 
large values. The VC flux was a particularly promising 
mechanism to export small-scale magnetic helicity because 
it could do so while allowing the system to retain much of 
the large-scale dynamo generated field, resulting in rapid 
growth to strong mean fields. Turbulent diff'usion of magnetic 
helicity is now the most promising escape from catastrophic 
a quenching - even for shearing systems. Simulation and 
theory suggest that this will be significant starting near 
ReM = lO'* (iMitra et al.l l20Tot Hubbard & Brandenburi 
1201 Oh . which, while astrophysically significant eludes numer- 
ical verification at present. However, diffusive fluxes must 
inevitably export scales of helicity at comparable fractional 
rates, reducing expected final field strengths below those that 
might have been hoped for with the VC flux. 
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